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The notion of singular reduction operators, i.e., of singular operators of nonclassical (conditional) 
symmetry, of partial differential equations in two independent variables is introduced. All possible 
' reductions of these equations to first-order ODEs are are exhaustively described. As examples, 

, properties of singular reduction operators of (1 -I- l)-dimensional evolution and wave equations are 

04 ' studied. It is shown how to favourably enhance the derivation of nonclassical symmetries for this 

^ , class by an in-depth prior study of the corresponding singular vector fields. 
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Distinctions in kind between Lie symmetries and nonclassical symmetries became apparent al- 
ready in the first presentation of nonclassical symmetries in [3| by the example of the (1 -|- 1)- 
dimensional linear heat equation and a particular class of operators. In contrast to classical Lie 
symmetries (see, e.g., [17]), the system of determining equations on the coefficients of nonclassi- 
cal symmetry operators of the heat equation was found to be nonlinear and less over deter mined, 
and the set of such operators does not possess the structure of an algebra or even a vector space. 
Another difference appears in the procedure of deriving the determining equations. Namely, 
' deriving systems of determining equations for nonclassical symmetries crucially depends on the 

interplay between the operators and the equations under consideration. Thus, for the linear heat 
i equation ut = Uxx the general form of nonclassical symmetry operators is Q = T{t,x,u)dt + 

I ^{t,x,u)dx + ri{t,x,u)du, where (r, ^) / (0,0), and there are two essentially different cases of 

. ' nonclassical symmetries: the regular case r 7^ and the singular case r = 0. The factorization up 

O 

00 : 1) r = 1 and 2) r = 0, C = 1. 

■ The problem of determining the nonclassical symmetries of the linear heat equation was 
completely solved in [10]. In the regular case r = 1, after partial integration of the corresponding 

I determining equations, we obtain ^ = g^{t,x) and ry = g'^{t,x)u -\- g^{t,x). The functions 

■ 9^ : 9^ and g^ satisfy a coupled nonlinear system of partial differential equations [3], which is 
linearized by a nonlocal transformation to a system of three uncoupled copies of the initial 
equation [9| IIOI [Mj . The underlying reason for this phenomenon lies in the interaction between 
the linearity and the evolution structure in the linear heat equation. Hence similar results can 
be obtained only for linear evolution equations [Tj [I8l [21] or related linearizable equations [13] . 

The singular case (r, ^) = (0,1) was not considered in [3]. In this case the system of de- 
termining equations for nonclassical symmetries consists of a single (1 -|- 2)-dimensional non- 
linear evolution equation for the unknown function 77 and, therefore, is not overdetermined. 
The determining equation is reduced by a nonlocal transformation to the initial equation with 
an additional implicit independent variable which can be assumed as a parameter |10] . The 
linearity of the heat equation is inessential here. Hence after the case of linear evolutions equa- 
tions [71 dH] this result was extended to general (1 -|- l)-dimensional evolution equations [25] . 
multi-dimensional evolution equations [19] and even systems of such equations [23]. Moreover, 
it was proved [H], that, e.g., in the (1 -|- l)-dimensional case there exists a one-to-one correspon- 
dence between one-parametric families of solutions of an evolution equation and its reduction 
operators with {t,£) = (0,1). 
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The above results raise a number of interesting questions, to wit: What are possible causes 
for the existence of singular cases for reduction operators? Is the conventional partition of sets of 
reduction operators with the conditions of vanishing and nonvanishing coefficients of operators 
universal or is it appropriate only for certain classes of differential equations, e.g., evolution 
equations? Can partitions of sets of reduction operators, different from the conventional one, be 
useful? Does there exist an algorithmic way of singling out singular cases for reduction operators 
before deriving determining equations? What properties of a partial differential equation and a 
subset of its reduction operators lead to a 'no-go' situation (i.e., a single determining equation 
equivalent, in a certain sense, to the initial equation)? What is the optimal way of obtaining 
the determining equation for nonclassical symmetries? The purpose of the present paper is to 
answer these and other related questions. 

Algorithms for deriving the determining equations for nonclassical symmetries were discussed, 
e.g., in [21 [6] but the focus of these works was quite different. 

The conditional invariance of a differential equation with respect to an operator is equivalent 
to any ansatz associated with this operator reducing the equation to a differential equation 
with one less independent variables [26]. That is why we use the shorter and more natural 
term "reduction operators" instead of "operators of conditional symmetry" or "operators of 
nonclassical symmetry" and say that an operator reduces a differential equation if the equation 
is reduced by the corresponding ansatz. The direct method of reduction with ansatzes of a 
special form was first explicitly applied in [3] to the Boussinesq equation although reductions 
by non-Lie ansatzes were already discussed, e.g., in [8j. A connection between the reduction 
by generalized ansatzes and compatibility with respect to higher-order constraints was found 
in [H]. 

To clarify the main ideas of the proposed framework of singular reduction operators, in this 
first presentation of the subject we consider only the case of a single partial differential equation 
in two independent and one dependent variables and a single reduction operator. We note, 
however, that more general cases can be included and will be the subject of forthcoming papers. 

Some of the main conclusions of the present paper are: 

• Singular cases of reduction operators of a partial differential equation are connected with 
the possibility of lowering the order of this equation on the manifolds determined by 
the corresponding invariant surface conditions in the appropriate jet space. Hence the 
first step of the procedure of finding nonclassical symmetries has to consist in studying 
singular modules of vector fields which lower the order of the equation. This step is 
entirely algorithmic, hence is especially suited to a direct implementation in symmetry- 
finding computer algebra programs. The structure of singular modules of vector fields has 
to be taken into account under splitting the set of reduction operators for factorization. 

• The weak singularity co-order of a reduction operator Q coincides with the essential order 
of the corresponding reduced equation and the number of essential parameters in the family 
of Q-invariant solutions. 

• If a single partial differential equation C in two independent variables admits a first co- 
order singular module S of vector fields then it necessarily possesses first co-order singular 
reduction operators belonging to S. The system of determining equations for such opera- 
tors consists of a single partial differential equation DE in three independent variables of 
the same order as C. The equation DE is reduced to £ by a nonlocal transformation. 

The paper is organized as follows: The main notions and statements on nonclassical symme- 
tries are presented in Section [2l Singular vector fields of differential functions and differential 
equations are defined and studied in Sections [3] and [H respectively. Singular reduction operators 
of (1 -|- l)-dimensional evolution and nonlinear wave equations are exhaustively investigated in 
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Sections [5] and [H It is shown that the conventional partition of sets of reduction operators is 
natural for evolution equations, in contrast to the case of nonlinear wave equations. A connec- 
tion between the singularity co-order of reduction operators and the number of parameters in 
the corresponding families of invariant solutions is established in Section [71 The final Section [8] 
is devoted to first co-order singular reduction operators of general partial differential equations 
in two independent and one dependent variables. 

2 Reduction operators of differential equations 

Following im [T2t l22l I26j . in this section we briefly collect the required notions and results 
on nonclassical (conditional) symmetries of differential equations. Also, we argue for the use 
of the name "reduction operators" instead of "nonclassical (conditional) symmetry operators" . 
In accordance with the aims of this paper we restrict our considerations to the case of two 
independent variables and a single reduction operator. 

The set of (first-order) differential operators (or vector fields) of the general form 

Q = e{x,u)di+vix,u)du, + (0,0), 

will be denoted by £}. In what follows, x denotes the pair of independent variables (xi,X2) and 
u is treated as the unknown function. The index i runs from 1 to 2, and we use the summation 
convention for repeated indices. Subscripts of functions denote differentiation with respect to 
the corresponding variables, di = d/dxi and du = d/du. Any function is considered as its 
zero-order derivative. All our considerations are carried out in the local setting. 

Two differential operators Q and Q are called equivalent if they differ by a multiplier which 
is a non-vanishing function of x and u: Q = XQ, where A = X{x,u), A 7^ 0. The equivalence of 
operators will be denoted by Q ~ Q. Factoring £2 with respect to this equivalence relation we 
arrive at £lf. Elements of £jf will be identified with their representatives in £3. 

The first-order differential function 

Q[u] := r]{x, u) - u)ui 

is called the characteristic of the operator Q. The characteristic PDF Q[u] =0 (also known 
as the invariant surface condition) has two functionally independent solutions C,{x,u) and 
Lo{x,u). Therefore, the general solution of this equation can be implicitly represented in the 
form -F(C, ^) = 0, where F is an arbitrary function. 

A differential function G = Q[z\ of the dependent variables z = (z^, . . . , z"*) which in turn 
are functions of a tuple of independent variables y = (yi, . . . , |/„) will be considered as a smooth 
function of y and derivatives of z with respect to y. The order r = ordO of the differential 
function Q equals the maximal order of derivatives involved in 0. More precisely, the differential 
function is defined as a function on a subset of the jet space J'^{y\z) [14] . 

The characteristic equations of equivalent operators have the same set of solutions. Con- 
versely, any family of two functionally independent functions of x and n is a complete set 
of integrals of the characteristic equation of a differential operator. Therefore, there exists a 
one-to-one correspondence between £lf and the set of families of two functionally independent 
functions of x and ti, which is factorized with respect to the corresponding equivalence relation. 
(Two families of the same number of functionally independent functions of the same arguments 
are considered equivalent if any function from one of the families is functionally dependent on 
functions from the other family.) 

Since {S},^^) 7^ (0,0) we can assume without loss of generality that C« / and 7^ and 
resolve the equation F = with respect io C,: Q = ^(uj). This implicit representation of the 
function u is called an ansatz corresponding to the operator Q. 
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Consider an rth order differential equation C of the form = for the unknown 

function u of two independent variables x = (xi,X2). Here L = L[u\ = L(x,U(r)) is a fixed 
differential function of order r and n(^) denotes the set of all the derivatives of the function 
u with respect to x of order not greater than r, including u as the derivative of order zero. 
Within the local approach the equation C is treated as an algebraic equation in the jet space 

_ jr^^i^-j q£ order r and is identified with the manifold of its solutions in J'': 

£ = {{x, n(j,)) G J*" I L(x, n(r)) = 0}. 

Denote the manifold defined by the set of all the differential consequences of the characteristic 
equation Q[u\ =0 in ,F by Q(r), i-e., 

Q(^) ={(x,U(^)) & ,r\D'^D'^Q[u] =0, a,/3GNU{0}, a + /?<r}, 

where Di = di + ^io+i,/?^^^^ and D2 = ^2 + Ua^p+idu^p are the operators of total differentiation 
with respect to the variables xi and 2:2, and the variable Uajs of the jet space corresponds to 

the derivative d'^'^^u/ 8x18x2- 

A precise and rigorous definition of nonclassical (or conditional) symmetry was first suggested 
in [n] (see also HZIES]). 

Definition 1. The differential equation C is called conditionally invariant with respect to 
the operator Q if the relation (5(j,)L(3;, M(j,))|^^g(,,)= holds, which is called the conditional 
invariance criterion. Then Q is called an operator of conditional symmetry (or Q-conditional 
symmetry, nonclassical symmetry, etc.) of the equation C. 

In Definition [T] the symbol (5(r) stands for the standard rth prolongation of the operator Q 

MM- 

Q{r)=Q+ J2 V""" ■=D'^D^Q[u]+eUa+i,p + fn^,0+i. 

0<a+/3^r 

The equation C is conditionally invariant with respect to Q if and only if the ansatz ^ = 99(0;) 
constructed with Q reduces C to an ordinary differential equation £: L[(p] = [26j. Namely, 
there exist differential functions A = X[ip] and L = L[ip] of an order not greater than r (i.e., 
functions of a; and derivatives of ip with respect to lj up to order r) such that L\u=^(uj) = 
The function A does not vanish and may depend on as a parameter, where the value 9 = 9{x,u) 
is functionally independent of C and lo. The differential function L is assumed to be of minimal 
order r which is possibly reached up to the equivalence generated by nonvanishing multipliers. 
Then the reduced equation C is of essential order f. 

This is why we will also call operators of conditional symmetry reduction operators of C. 

Another treatment of conditional invariance is that the system C Ci Q(j.) is compatible in the 
sense of not involving any nontrivial differential consequences \15\ [16] . 

The property of conditional invariance is compatible with the equivalence relation on Q 

Lemma 1. // the equation C is conditionally invariant with respect to the operator Q then it is 
conditionally invariant with respect to any operator which is equivalent to Q. 

The set of reduction operators of the equation £ is a subset of and so will be denoted by 
£l(£). In view of LemmalU Q G £2(£) and Q ^ Q imply Q £ £2(£), i.e., 0.{C) is closed under 
the equivalence relation on 0. Therefore, the factorization of with respect to this equivalence 
relation can be naturally restricted to 0(£), resulting in the subset 0f (£) of 0f . As in the whole 
set 0f, we identify elements of 0.f{C) with their representatives in £}.{C). In this approach the 
problem of completely describing all reduction operators for C is equivalent to finding 0f(£). 

The conditional invariance criterion admits the following useful reformulation [26] . 
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Lemma 2. Given a differential equation L: L[it] = of order r and differential functions L[u] 
and \[u\ ^ of an order not greater than r such that L\q,^^ = ^LlGr^.)) o,n operator Q is a 



The classification of reduction operators can be considerably enhanced and simplified by 
considering Lie symmetry and equivalence transformations of (classes of) equations. 

Lemma 3. Any point transformation of x and u induces a one-to-one mapping of O. into itself. 
Namely, the transformation g: Xi = X^{x^u), u = U{x,u) generates the mapping : £} — > H 
such that the operator Q is mapped to the operator g^Q = Cdxi+V^u, where ^{x^ u) = QX^{x, u), 
fj{x,ii) = QU{x,u). If Q' ^ Q then g^,Q' ~ g^^Q. Therefore, the corresponding factorized 
mapping gf.£}.f^ £2f also is well defined and bijective. 

Definition 2 ([20]). Differential operators Q and Q are called equivalent with respect to a 
group G of point transformations if there exists g € G for which the operators Q and g^Q are 
equivalent. We denote this equivalence by Q ~ Q mod G. 

Lemma 4. Given any point transformation g of an equation C to an equation C, maps Q.{C) 
to Q.{C) bijectively. The same is true for the factorized mapping gf from 0-{{C) to £)f(£). 

Corollary 1. Let G be the point symmetry group of an equation C. Then the equivalence of 
operators with respect to the group G generates equivalence relations in il{C) and in l}f(/2). 

Consider the class £1^ of equations Cg: L(j;, 6*) = parameterized with the parameter- 
functions = 9{x,U(^j.))- Here L is a fixed function of x, M(^) and 6. The symbol 9 denotes 
the tuple of arbitrary (parametric) differential functions d{x, M(r.)) = {9^{x, n(r)), . . . , d'^{x, «(,■))) 
running through the set S of solutions of the system 5(a;, ■U(r), 0(q)(x, U(r))) = 0. This system 
consists of differential equations on 9, where x and play the role of independent variables 
and 0(g) stands for the set of all the derivatives of of order not greater than q. In what follows 
we call the functions arbitrary elements. Denote the point transformation group preserving 
the form of the equations from jC\s hy G^ . 

Let P denote the set of the pairs consisting of an equation Cq from £1^ and an operator Q 
from 0.{Ce). In view of Lemma HJ the action of transformations from the equivalence group G~ 
on C\s and {'lI{Cq) | G 5} together with the pure equivalence relation of differential operators 
naturally generates an equivalence relation on P. 

Definition 3. Let 9,9' eS,Q e Q{Cg), Q' G £l{Ce'). The pairs {Ce,Q) and {Ce',Q') are called 
G^ -equivalent if there exists g G G~ such that g transforms the equation Cg to the equation jCg^, 
and Q' ~ g^Q. 

The classification of reduction operators with respect to G~ will be understood as the clas- 
sification in P with respect to this equivalence relation, a problem which can be investigated 
similar to the usual group classification in classes of differential equations. Namely, we construct 
firstly the reduction operators that are defined for all values of 9. Then we classify, with respect 
to the values of 9 for which the equation Cq admits additional reduction operators. 

3 Singular vector fields of differential functions 

Consider a vector field Q = Ci^^ u)di-\-r]{x, u)du with {£} 7^ (0, 0), defined in the space (x, u), 
and a differential function L = L[u\ of order ordL = r (i.e., a smooth function of x = (xi, 2:2) 
and derivatives of u of orders up to r). 




holds, where r 



ordL ^ r 
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Definition 4. The vector field Q is called singular for the differential function L if there exists 
a differential function L = L\u] of an order less than r such that L\qi^^.^ = L|qj^,j. Otherwise Q 
is called a regular vector field for the differential function L. If the minimal order of differential 
functions whose restrictions on Q(^) coincide with L\q^^^^ equals k (k < r) then the vector field Q 
is said to be of singularity co- order k for the differential function L. The vector field Q is called 
ultra- singular for the differential function L if L\Q(r) = ^■ 

For convenience, the singularity co-order of ultra-singular vector fields and the order of iden- 
tically vanishing differential functions are defined to equal —1. Regular vector fields for the 
differential function L are defined to have singularity co-order r = ordL. The singularity co- 
order of a vector field Q for a differential function L will be denoted by sco^ Q. 

If Q is a singular vector field for L then any vector field equivalent to Q is singular for L with 
the same co-order of singularity. 

A function L satisfying the conditions of Definition |4] can be constructively found. Namely, 
without loss of generality we can suppose that the coefficient of 82 in Q is nonzero. Then any 
derivative of u of order not greater than r can be expressed, on the manifold Q(r)) via derivatives 
of u with respect to xi only. For example, for the first- and second-order derivatives we have 

U2 = 'n- iui, 

U12 = m - iiui + f]uUi - ^uuj - iuu, (1) 

U22 =m- i2Ul + {fju - iuUl){fj - |ui) - 1(571 - ^lUi + fjuUl - iuul - iuil), 

where ^ = ^^/^^ and f/ = rj/^^ . After substituting the expressions for the derivatives into L, we 
obtain a differential function L depending only on x, u and derivatives of u with respect to xi. 
We will call L a differential function associated with L on the manifold Q(r)- The vector field Q 
is singular for the differential function L if and only if the order of L is less than r. The co-order 
of singularity of Q equals the order of L. The vector field Q is ultra-singular if and only if L = 0. 
Therefore, testing that a vector field is singular for a differential function with two independent 
variables is realized in an entirely algorithmic procedure and can be easily included in existing 
programs for symbolic calculations of symmetries. 

Consider the two-dimensional module {Q^ = 6^Q^} of vector fields over the ring of 
smooth functions of {x,u) generated by the vector fields = ^^^{x,u)dj + r]^{x,u)du, where 
rank(^*^, (^*^, r/*) = 2. In the remainder of this section the parameter tuple 9 = {9^,6'^) runs 
through the set of pairs of smooth functions depending on (x,n), and i and j run from 1 to 2. 

Definition 5. The module {Q^} is called singular for the differential function L if for any 9 
with (^*^*^, ^*^*^) 7^ (0,0) the vector field is singular for L. The singularity co-order of the 
module {Q^} coincides with the maximum of the singularity co-orders of its elements. 

By a point transformation, one of the basis vector fields, e.g. Q^, can be reduced to du (trans- 
forming L simultaneously with and Q^.) Then {^^^iS}"^) 7^ (0,0), and up to permutation of 
independent variables we can assume 7^ and, therefore, set r/^ = and S}"^ = 1 by a change 
of basis. Any vector field from the module {Q^} with a nonzero value of 9^ is equivalent to the 
vector field + CQ^) where C, = 9"^ /9^ . All the other vector fields from {Q^} (which have 0^ = 
and, therefore, are equivalent to du) can be neglected since each of them leads to the equation 
9'^{x,u) = which completely determines u and therefore, does not give an ansatz for u. 

This justifies why, up to point transformations, it suffices to study only singular sets of vector 
fields of the form {Q'^ = ^di + ^2 + C^w}; with ^ a fixed smooth function of (x, u) and C, running 
through all such functions. The latter form of singular sets of vector fields will be called reduced. 

Further simplification depends on whether the module is closed under the Lie bracket. In 
case it is, it can be assumed to be generated by two commuting vector fields which can be 
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simultaneously reduced by a point transformation to shift operators, e.g., = 82 and = d^- 
Thus in the reduced form ^ can be put to 0. If the module is not closed under the Lie bracket, 
we have 7^ in the reduced form. After the point transformation Xj = Xj and = ^ and a 
change of basis, we obtain the basis = udi + 82 and = du. Hence: 

Proposition 1. In any two-dimensional module of vector fields in the space of three variables 
{xi,X2,u), any basis vector fields and can be locally reduced, by point transformations, to 
the form = 82 (resp. = udi + 82) and = 8u if the module is closed (resp. not closed) 
with respect to the Lie bracket of vector fields. 

Theorem 1. A differential function L with one dependent and two independent variables pos- 
sesses a kth co-order singular two-dimensional module of vector fields if and only if it can be 
represented, up to point transformations, in the form 

L = L{x,Qr,k), (2) 

where Qj.,k = {^a = {^Di + D2)'^^u, ai ^ k,ai + 02 ^ r^, ^ S {0, u}, and L^^ 7^ for some 
uja with ai = k. 

Proof. Suppose that a differential function L possesses a A:th co-order singular two-dimensional 
module of vector fields {Q^ = 9^Q^}. By a point transformation and a change of basis, we 
represent the basis elements in the reduced form = -|- 82 and = 8u, where ^ G {0, «}, 
and choose the subset {Q^ = -\- 82 + Cdu} in {Q^}, where C runs through the set of smooth 
functions of (x, u). The initial differential function also will be changed by these transformations 
but throughout we will use the old notations for all new values. 

We fix an arbitrary point = {x^,u1^^^) G and consider the vector fields from {Q''} for 
which G '^^^^ condition implies that the values of the derivatives of C with respect 

to only xi and X2 in the point (x^,u^) are expressed via u^^-j and values of derivatives of C in 
{x^,u^), containing differentiation with respect to u. The latter values are not constrained. 

We introduce the new coordinates {xj, Wq, = D"^ {(.Di + D2)'^'^u, \a\ ^ r} in instead of the 
standard ones {xj,«Q, |a| ^ r}. This is a valid change of coordinates since the Jacobian matrix 
{8uja/8uai) is nondegenerate. Indeed, it is a triangular matrix with all diagonal entries equal to 1 
if the following order of multi-indices is implemented: a < (3 -.^^ \a\ < \(3\ V (|a| = \(3\ Aa2 < (i2)- 
Note that u;„ = D'^^^Di + D2Y^u = D'^^iQ^y^u on Q[ 

Consider the differential function L obtained from L by the above procedure of excluding, 
on the manifold Q/^jy the derivatives of u involving differentiations with respect to X2 (see ([1])). 
Since Q'^ is a kth co-order singular vector field for L, the function L does not depend on the 
derivatives M(k,o)) k = 1, . . . , r. We use this condition step-by-step, starting from the greatest 
value of K and re- writing the derivatives in the new coordinates of J*" and in terms of L. 

Thus, in the new coordinates the equation ^^^^^^(z'^) = has the form ^^^^^^(z'^) = 0. This 

completes the first step. Then in the second step the equation (z*^) = implies that 

We split with respect to the value Cu{x^,u^) since it is unconstrained. As a result, we obtain 
the equations L^^^^^^^^^{z^) = and L^^^_^^^^{z^) = 0. 

Iterating this procedure, before the yuth step, ^ G {1, . . . ,r — A;}, we derive the equations 
^^(r-M',^) = 0, /x' = 0, . . . , /i - 2, 1/ = 0, . . . , Then the equation {z^) = implies 

that 
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The values d'^'^^C{x^ ,u^), = 0, — 1, are unconstrained. Splitting with respect to them, 
which is equivalent to splitting with respect to {du{Q'^Yu) ^)=(2,o ^oy = 0, . . . , — 1, gives 
the equations ^a;(,_^+i,,) {z^) = 0, = 0, . . . , - 1. 

Finally, after the (r — A;)th step we derive the system L^^^_^, (z^) = 0, fi' = 0, . . . ,r — k + 1, 
u = 0, . . . , n' , which implies condition ([2]). 

Conversely, let an rth order differential function L be of the form ([2]) (after a point 
transformation). For an arbitrary smooth function ( = ({x,u) we consider the vector field 

= £^di + ^2 + C9u and the differential function L = L{x, 0,r,k) where 

^r,k = (a;„ = D^^(Q^)°2n,Qi ^ A:, ai + aa < r) . 
Then ord L = k and 



i.e., {Q^ = + CQ^}; where = ^di + 82, = du and C runs through the set of smooth 
functions of {x,u), is a kth co-order singular set for the differential function L in the new 
variables. We complete the set by the vector fields equivalent to its elements or du and return to 
the old variables. As a result, for the differential function L we construct a kth co-order singular 
two-dimensional module of vector fields {Q^ = 6^Q^}. □ 

Corollary 2. A differential function with one dependent and two independent variables admits 
a kth co-order singular two-dimensional module generated by commuting vector fields if and only 
if it can be reduced by a point transformation of the variables to a differential function in which 
all differentiations with respect to one of the independent variables are of order < k. 

Corollary 3. Any differential function with one dependent and two independent variables (not 
identically vanishing) admits no ultra- singular two-dimensional module of singular vector fields. 

Note 1. It is obvious that a singular module may contain vector fields whose singularity co- 
orders are less than the singularity co-order of the whole module. Suppose that {Q^ = (,di -\- 
d2 + Cdu} is a singular set of vector fields for a differential function L, and its singularity co-order 
equals k. Then the values of C for which sco^ < k are solutions of the equation 

r—k 

Y,L^^,,^,{xA,k){du{Q^ru) =Q, 
where 0^,/, = (L»^i(g«)°2 

u, ai ^ k,ai + a2 ^ r) and L is defined in Theorem[TJ In other words, 
the regular values of C, associated with the vector fields of the maximal singularity co-order k in 
{Q''} satisfy the inequality 

r—k 

Y,L^^^^^{x,^r,k){du{Q^ru) ^Q. 

4 Singular vector fields of differential equations 

We will say that a vector field Q is (strongly) singular for a differential equation C if it is singular 
for the differential function L[u] which is the left hand side of the canonical representation 
L[u] =0 of the equation C Usually we will omit the attribute "strongly". 

Since left hand sides of differential equations are defined up to multipliers which are nonvan- 
ishing differential functions, the conditions from Definition S] can be weakened when considering 
differential equations. 
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Definition 6. A vector field Q is called weakly singular for the differential equation C: L[n] = 
if there exist a differential function L = L[u\ of an order less than r and a nonvanishing differential 
function A = A[n] of an order not greater than r such that L\q^^^ = \L\ q^^^ . Otherwise Q is called 
a weakly regular vector field for the differential equation C If the minimal order of differential 
functions whose restrictions on Q(j,) coincide, up to nonvanishing functional multipliers, with 
L\Q(r) equal to k {k < r) then the vector field Q is said to be weakly singular of co-order k for 
the differential equation C 

The notions of ultra-singularity in the weak and the strong sense coincide. Analogous to the 
case of strong regularity, weakly regular vector fields for the differential equation C are defined 
to have weak singularity co-order r = ordL. The weak singularity co-order of a vector field Q 
for an equation C will be denoted by wsco£ Q. 

Note that strong singularity implies weak singularity and hence weak regularity implies strong 
regularity. It is obvious that the weak singularity co-order is never greater and may be less than 
the strong singularity co-order. In particular, strongly regular vector fields may be singular in 
the weak sense. For example, the equation uut = e""""" {ux + u) possesses the singular vector field 
dt whose strong and weak singularity co-order equal 2 and 1, respectively. The same vector field 
dt is strongly regular and is of weak singularity co-order 1 for the equation ut = e^'^''{ux + u). 

If Q is a weakly singular vector field for C then any vector field equivalent to Q is weakly 
singular for C with the same co-order of weak singularity. 

Weakly singular vector fields are related to characteristic directions (cf. |14j concerning char- 
acteristic directions and characteristic hypersurfaces): Given a vector field Q = ^^{x,u)di + 
r]{x,u)du weakly singular for a differential equation £, in each point of the manifold C the 
vector (^^,^^) is orthogonal to a characteristic direction of the equation £ in this point. 

Let L be a differential function associated with L on the manifold Q(r), namely, obtained 
from L via excluding those derivatives of u which contain differentiations with respect to X2 
in view of equations defining Q(r)- Suppose additionally that L is of maximal rank in the 
derivative u of the highest order k appearing in this differential function, i.e., o) 7^ on 
the solution manifold of the equation L = 0. Then the weak singularity co-order of Q for the 
equation C: L = equals the order k of L and, therefore, the strong singularity co-order of Q. 
Hence in this case testing that a vector field is weakly singular for a partial differential equation 
with two independent variables can be implemented via an entirely algorithmic procedure. 

Theorem 2. An rth order differential equation C: L[u] = of maximal rank with one dependent 
and two independent variables possesses a kth co-order weakly singular two-dimensional module 
of vector fields if and only if L can be represented, up to point transformations, in the form 

L = A[u]L{x,nr,k), (3) 

where A is a nonvanishing differential function of order not greater than r, Qj.,k = (<^o = 
L*"^ {S,Di + D2)'^^u, ai ^ /c, ai + a2 ^ ; C £ {0, u}, and L^^ ^ for some uJa with ai = k. 

Proof. We will freely use the notations and definitions from the proof of Theorem [TJ 

Suppose at first that a differential equation C: L[u] = is of maximal rank and admits a kth 
co-order weakly singular two-dimensional module of vector fields. Up to point transformations 
and changes of module basis, we may consider only a set {Q^ = ^di + ^2 + (du} of singular 
vector fields in reduced form. 

We fix an arbitrary point = (x^,u9 n) G £ C J*" and choose the vector fields from {Q''} 
for which z G 2(r)' This condition implies that the values of derivatives of C with respect to 
only xi and X2 in the point {x^,u^) are expressed via li^^^ and values of derivatives of C, in 
{x^,u^), containing differentiation with respect to u. The latter values are not constrained. 
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The differential function L is obtained from L by excluding, on the manifold Q^^^ , derivatives 
of u involving differentiations with respect to X2 (see ([IJ). kth co-order weak singularity of 
for L leads to mi^o) = 0, k = A; + 1, . . . , r. We use this condition step-by-step as in the 
proof of Theorem [H starting from the greatest value of k and re- writing the derivatives in the 
new coordinates {xi^LOa = + D2)"^u, \a\ ^ r} of J*" and in terms of L. Therefore, 

=0' /"' = 0,...,r-/c + l, z. = 0,...,/x', 

which is satisfied for any £ C Applying the Hadamard lemma to each of these equations and 
then simultaneously integrating them, we obtain ([3]) (cf. the proof of Theorem 1 in [26j). 

Conversely, let an rth order differential function L be of the form ([3]) (after a point 
transformation). For an arbitrary smooth function C = ("(x,n) we consider the vector field 

= + ^2 + C^u and the differential function L = L{x, 0.r,k)-, where 

^r,k = {uj^ = (Q'^T'u, ai ^ fc, ai + a2 ^ r) . 

ThenordL = k and L|„c = A^l^c > i-e-, IQ^ = Q^ + CQ^I, where = ^81 + 82, = 8u and 
^ runs through the set of smooth functions of {x,u), is a kth co-order weakly singular set for the 
differential equation C in the new variables. We complete the set by the vector fields equivalent 
to its elements or 8u and return to the old variables, thereby constructing a kth co-order weakly 
singular two-dimensional module of vector fields {Q^ = O^Q^} for the differential equation C. □ 

Corollary 4. A differential equation C: L[u] = of maximal rank with one dependent and two 
independent variables possesses a kth co- order weakly singular two-dimensional module of vector 
fields if and only if this module is kth co-order strongly singular for C (possibly in a representation 
differing from L[u] = in multiplication by a nonvanishing differential function of u). 

Definition 7. A vector field Q is called a singular reduction operator of a differential equation C 
if Q is both a reduction operator of C and a weakly singular vector field of C 

5 Example: evolution equations 

In this section we investigate singular reduction operators of (1 + l)-dimensional evolution equa- 
tions of the form 

Ut = H{t,X,Ui^r,x)), (4) 

where r > 1, uq := u, Uk = 8^u/dx^, ^(r,x) = {uq^ui, . . . ^Ur) and Hu^ 7^ 0. (We revert to 
the notation t and x for xi and X2, respectively, and change the notations of the corresponding 
derivatives.) Evolution equations are quite specific from the point of view of singular vector 
fields and singular reduction operators. 

Proposition 2. A vector field Q = T{t, x,u)8t + ^{t,x,u)8x -\- r]{t, x,u)8u is singular for the 
differential function L = ut — H{t, x, 0/ order r > 1 if and only if t = 0. The co-order of 

singularity of any singular vector field for any such differential function equals 1. 

Proof. Suppose that r 7^ 0. Excluding the derivative ut from L according to the equation 
Ut = rj/r — S,Ux/t results in a differential function L = rijr — ^Uxjr — H(t, x,U(^r,x))- Since 
ordL = r = ordL, the vector field Q is not singular in this case. Therefore, for the vector field 
Q to be singular, the coefficient r has to vanish. 

If r = and therefore 7^ 0, all the derivatives n^t, k = 1, . . . ,r, can be expressed, on the 
manifold Q(r) via t, x and u: = {8x + C,du)^~^C,, k = 1, . . . ,r, where C, = r]/£_. Using these 
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expressions for excluding the derivatives n^, k = l,...,r from L, we obtain the differential 
function 

L = ut- H{t, x,u), H := H{t, x, u, (, Cx + CC«, • • • , {dx + Cduf-K), 

whose order equals 1. Hence the vector field Q is singular for the differential function L, and its 
singularity co-order equals 1. □ 

Corollary 5. For any (1 + \)- dimensional evolution equation, the corresponding differential 
function possesses exactly one set of singular vector fields in the reduced form, namely, S = 
{dx + C{x,u)du}- The singularity co-order of S equals 1. 

It is obvious that under the condition Hu^ 7^ a vector field is singular for the differen- 
tial function ut — H{t^x,u^j.„^~^) if and only if it is weakly singular for the differential equation 
ut = Hit^x^Uf^^;^-^). Hence we do not distinguish between strong and weak singularity (cf. 
Corollary H]). 

The vector fields 82 and du generating the singular module associated with 5 commute and 
the differential function L contains only first order differentiation with respect to t (namely, in 
the form of the derivative ut). This perfectly agrees with Corollary [2l 

We fix an arbitrary equation C of the form @ and denote by Qo{C) the set of reduction 
operators of £, belonging to S. For the equation C and Q G Qq{C)^ the conditional invariance 
criterion implies only the single rth order equation 

Ct + CuH = Hx + A, H := H{t, X, u, (, Cx + CCu, ■ ■ ■ , {dx + CduY-\), 

with respect to the single unknown function C with three independent variables t, x and u, 
which we will denote by DEo(>C). In other words, the system of determining equations in this 
case consists of the single equation DEo(>C) and, therefore, is not over determined. DEq(C) is the 
compatibility condition of the equations = C and C. 

Theorem 3. Up to the equivalences of operators and solution families, for any equation of 
form dH) there exists a one-to-one correspondence between one-parametric families of its solu- 
tions and reduction operators with zero coefficients of dt- Namely, each operator of this kind 
corresponds to the family of solutions which are invariant with respect to this operator. The prob- 
lems of the construction of all one-parametric solution families of equation ([4]) and the exhaustive 
description of its reduction operators with zero coefficients of dt are completely equivalent. 

Proof. Let C be an equation from class ^ and Q = dx -\- Qdu G Qq{C), i.e., the coefficient 
C = C,{t,x,u) satisfies the equation DEo(>C). An ansatz constructed with Q has the form u = 
f{t,x,(piuj)), where / = f{t,x,ip) is a given function, /(^ 7^ 0, (/J = ip{uj) is the new unknown 
function and uj = t is the invariant independent variable. This ansatz reduces £ to a first- 
order ordinary differential equation C in ip, solvable with respect to (p' . The general solution 
of the reduced equation C can be represented in the form = 99(0;, x), where ip^ ^ and x is 
an arbitrary constant. Substituting this solution into the ansatz results in the one-parametric 
family of solutions u = f{t, x, x) of C with / = /(t, x, (p{t, x)). Expressing the parameter x 
from the equality u = f{t,x,x), we obtain that x = ^{t,x,u), where 7^ 0. Then C = = 
—^x/^u for any u £ T, i.e., for any admissible value of (t, x, x). This implies that C = —(^x/^u 
for any admissible value of {t,x,u). 

Conversely, suppose that = {u = f{t,x, x)} is a one-parametric family of solutions of C. 
The derivative f^ is nonzero since the parameter x is essential. We express x from the equality 
u = f{t,x,x): X = ^(t,x,u) for some function $ = ^{t,x,u) with 7^ 0. Consider the 
operator Q = dx + (du, where the coefficient ( = C{t,x,u) is defined hy ( = —^x/^u- Q[u] = 
for any u £ T. The ansatz u = f{t,x,ip{uj)), where uj = t, associated with Q, reduces £ to the 
equation ip^ = 0. Therefore [26], Q G Qo{C) and hence the function ( satisfies DEq(jC). □ 
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Corollary 6. The nonlinear (1 + 2) -dimensional evolution equation DEo(£) is reduced by the 
composition of the nonlocal substitution = —^x/^u, where ^ is a function of {t, x, u), and the 
hodograph transformation 

the new independent variables: t = t, x = x, x = 
the new dependent variable: u = u 

to the initial equation C in the function u = ii{t,x,>c) with x playing the role of a parameter. 

Note 2. One-parametric families u = f{t,x,>c) and u = f{t,x,>c) are defined to be equivalent 
if they consist of the same functions and differ only by parameterizations, i.e., if there exists a 
function ( = C,{>c) such that C>c / and f{t,x,({^)) = f{tiX,>c). Equivalent one-parametric 
families of solutions are associated with the same operator from Qo{C) and have to be identified. 

Note 3. The triviality of the above ansatz and the reduced equation results from the above 
special representation for the solutions of the determining equation. Under this approach diffi- 
culties in the construction of ansatzes and the integration of the reduced equations are replaced 
by difficulties in obtaining the representation for the coefficients of the reduction operators. 

The above consideration shows that for any evolution equation C the conventional partition 
of the set £}.{C) of its reduction operators with the conditions r 7^ and r = is natural 
since it coincides with the partition of into the singular and regular subsets. This is 

a specific property of evolution equations which does not hold for general partial differential 
equations in two independent variables. After factorizing the subsets of £l{C) with respect to 
the usual equivalence relation of reduction operators, we obtain two different cases of inequivalent 
reduction operators (the regular case r = 1 and the singular case r = and ^ = 1), which have 
to be studied separately. 

Singular reduction operators of C are described in a unified 'no-go' way. All singular reduction 
operators of C have the same singularity co-order equal to 1 and hence reduce C to first order 
ordinary differential equations. The coincidence of the singularity co-orders guarantees the 
existence of a bijection between the set of singular reduction operators of £ and the set of one- 
parametric families of its solutions (up to the natural equivalence relations in these sets). As 
a result, in the case r = and = 1 the determining equation for a single coefficient of du is 
reduced, with no additional assumptions and conditions, to the initial equation Chy a nonlocal 
transformation (cf. Corollary [6l) . 

The regular case r = 1 is more complicated than the singular one. It essentially depends 
on the structure of the equation including the order, the kind of nonlinearities, etc. Up to now 
there are no exhaustive results on regular reduction operators even for second-order evolution 
equations. Only certain subclasses of such equations were investigated. See, e.g., [U [Sj [211 [22] for 
the complete classifications of regular reduction operators for some subclasses of second-order 
evolution equations parameterized by functions of single arguments. For example, even for the 
class of nonlinear diffusion equations of the general form ut = {f{u)ux)x (a classical example of 
solving a group classification problem for partial differential equations [T7] ) , the set of values of 
the parameter-function / which correspond to equations possessing non-Lie regular reduction 
operators has not yet been found. Most evolution equations have no regular reduction operators. 
A simple example is 

Ut = Uxx + we"" + xe^"" + ie^"=" + e"^"" + e^"*" . 

Some evolution equations (the linear ones [lOlllI], Burgers' equation [13], etc.) possess so many 
regular reduction operators that 'no-go' statements like those for singular reduction operators 
are true for them, but the nature of this 'no-go' differs from the 'no-go' of the singular case and 
is related to the property of linearity or linearizability of the corresponding evolution equations. 
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6 Example: nonlinear wave equations 



The next example which we study in detail within the framework of singular reduction operators 
is given by the class of nonlinear wave equations (in the characteristic, or light-cone, variables) 
of the general form 

ui2 = F{u). (5) 

Here F is an arbitrary smooth function of u. This class essentially differs from the class of 
evolution equations within the framework of singular vector fields. The main differences are 
that each differential function corresponding to an equation from class ([5]) has two singular 
sets of vector fields and these sets contain vector fields of lower singularity co-orders than the 
singularity co-orders of the whole sets. Thus, for any F the vector field Q = u)di+r]{x, u)du 
is singular for the corresponding differential function L = ui2 — F{u) if and only if ^^^^ = 0. 
Moreover, it is obvious that there are no differences between strong and weak singularity of 
vector fields for equations from the class dS]). Indeed, suppose that 7^ 0. Excluding the 
derivatives U2 and U12 from L according to ([T]), we obtain a differential function L with the 
coefficient —£} H'^ of un. We have ordL < 2 if and only if = 0. 

Therefore, for any F the differential function L = U12 — F{u) possesses exactly two sets of 
singular vector fields in the reduced form, S = {82 + Ci^^ u)du} and S* = {di + C*{^^ u)du}- The 
vector fields equivalent to du are not suitable as reduction operators. Any singular vector field 
of L is equivalent to one of the above fields. Moreover, each equation of the form ^ admits 
the discrete symmetry transformation permuting the variables xi and X2- This transformation 
generates a one-to-one mapping between S and 5* (cf. Corollary [1]) . Hence it suffices, up to 
equivalence of vector fields (and permutation of xi and X2), to investigate only singular reduction 
operators from the set S. 

For an equation £ from class (jS]) and an operator Q = 82 + C^u the conditional invariance 
criterion takes the form 

(C12 + CluU2 + C,2uUl + CuuUlU2 + C«'"12) | £nS(2) = (Fu- 

The intersection C n Q(2) is singled out from by the equations U2 = C) Ci + CuUi = F and 
U12 = F. Our further considerations therefore depend on the values of Cm and Fu- We analyze 
all the possible cases. 

Let Cm = aiid F^ = 0. Then Q is an ultra-singular vector field for the differential function L. 
The third equation defining C Ci Q(2) takes the form (^i = F and contains no derivatives of u. 
It should be assumed as a condition with respect to C and hence the conditional invariance 
criterion is identically satisfied in this case. An ansatz constructed with the operator Q is 
u = (p{uj) + / Cdx2, where u; = xi. It reduces equation ([5]) to an identity. This is explained by 
the ultra-singularity of the reduction operator Q. 

If Cm = and Fu ^ 0, the singularity co-order of Q for the differential function L equals 0. 
The third equation defining C D Q(2) again takes the form Ci = but now can be solved with 
respect to u: u = F{(i), where F is the inverse to F. Then the conditional invariance criterion 
is equivalent to the equation C12 = C-^m(-^(Ci)) with respect to C- The ansatz constructed with 
the operator Q reduces equation ([5]) to the algebraic equation F{if + jC,dx2) = Ci for the 
function if in agreement with the singularity co-order of Q. Indeed, inverting F, we obtain 
the equality (p = F{(i) — J C,dx2 whose right-hand side does not depend on X2 in view of the 
equation on C- Conversely, let us fix a solution u = f{x) of equation ([5]) and set C = /2- Then 
C12 = C-^m(-^(Ci))> i-e-, in view of the conditional invariance criterion Q = d2 + Cdu is a reduction 
operator of equation ([5|), and Cm = 0. The solution u = f{x) is invariant with respect to Q. The 
above results can be summed up as follows: 
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Theorem 4. For any equation from class ([5]) with Fy_ ^ there exists a one-to-one corre- 
spondence between its solutions and reduction operators of the form Q = 82 + C{x)du (resp. 
Q* = di -\- C,*{x)du)- Namely, each operator of this kind is of singularity co-order and corre- 
sponds to the solution which is invariant with respect to this operator. The problems of solving 
an equation from class ([5]) with Fu ^ and the exhaustive description of its reduction operators 
of the above form are completely equivalent. 

Corollary 7. Any solution u = f{x) of equation ([5]) with F^ ^ is invariant with respect to two 
reduction operators Q = 82 + C{x)du and Q* = di + C*{x)du of equation ([5]), having singularity 
co-order 0. Here C = /2 and C* = /i- property of possessing the same invariant solution of 
equation ([5]) establishes a canonical bijection Q ^ Q* between the sets of reduction operators of 
singularity co-order 0. The adjoint values of C, and Q* are connected by the formulas 

^* _ Cii ^ _ C22 



F„(F(Ci)) Fu{F{Q)) 

The regular values of C, for which the singularity co-order of Q coincides with the singularity 
co-order of the whole family S (and equals 1) satisfy the condition Qu 7^ 0. The third equation 
defining C n Q(2) then provides the following expression for m: 

ui = — =: C ■ 

The conditional invariance criterion implies only the single equation 

C12 + CClu + {C2u + CCuu)^^ +CuF = (Fu (6) 

with respect to the single function i.e., in this case the system of determining equations 
consists of the single equation ^ and, therefore, is not over deter mined. 
Equation ^ can be rewritten in the form of the compatibility condition 

Ci + CCu = C2 + CCu = F 

of the equations ui = C*, U2 = C ^12 = F. It is obvious that C* 7^ 0. Due to symmetry with 
respect to the permutation of xi and X2, we obtain the following statement. 

Proposition 3. For any equation from class ([5]), there exists a canonical bijection Q <^ Q* 
between sets of its singular reduction operators of the forms Q = 82 + Cix,u)du and Q* = 
9i + C*(2^i u)du, where (^^ ^ and (* 7^ 0. This bijection is given by the formulas 

A solution of equation ([5]) is invariant with respect to the operator Q if and only if it is invariant 
with respect to the operator Q* . 

Theorem 5. Up to the equivalence of solution families, for any equation from class ([5]) with 
Fu ^ there exists a one-to-one correspondence between one-parametric families of its solutions 
and reduction operators of the form Q = d2+C{x, u)du, where Cu 7^ (resp. Q* = di+C*{x,u)du, 
where C* 7^ 0). Namely, any such operator corresponds to the family of solutions which are 
invariant with respect to this operator. The problems of the construction of all one-parametric 
solution families of an equation from class ([5]) with F^ ^ and the exhaustive description of its 
reduction operators of the above form are completely equivalent. 
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Proof. In view of Proposition [3l it is sufficient to consider only operators with zero coefficient of 
di . Although the proof is similar to the proof of the analogous statement for evolution equations 
it differs from it in essential details and will therefore be presented completely. 

An ansatz constructed with the operator Q = 82 + C,{x,u)du has the form u = f{x,ip{uj)), 
where / = f{x,ip) is a given function, 7^ 0, 93 = f{uj) is the new unknown function and 
(J = xi is the invariant independent variable. Here 7^ implies 7^ 0. Hence this ansatz 
reduces equation ([5]) to a first-order ordinary differential equation C in ip, which is solvable 
with respect to ip'. The general solution of the reduced equation C essentially depends on an 
arbitrary constant x: tp = ip{io,>c), where ip;^ 7^ 0. Substituting the general solution into the 
ansatz gives the one-parametric family T of solutions u = f{x, x) of ([5]) with / = f{x, ip{xi, x)). 

Conversely, suppose that F„ / and = {u = is a one-parametric family of 

solutions of ([5]). The derivative is nonzero since the parameter x is essential. Therefore, 
fi2}< = f}<Fuif) 7^ 0. We express x from the equality u = /(x,x): x = <I>(x,u) for some 
function $ = ^{x,u) with 7^ 0. Consider the operator Q = 82 + (du, where the coefficient 
C = Ci^ju) is defined by the formula C = QH = for any u £ T. The ansatz 

u = f{x,ip{uj)), where lj = xi, associated with Q, reduces ([5]) to the equation Pui = since 
/2x 7^ 0. Therefore [26], Q is a reduction operator of equation ^ and hence the function ^ 
satisfies equation ([6]). Moreover, we have Cu 7^ since otherwise the operator Q would reduce ([5]) 
to an algebraic equation with respect to (p. □ 

Corollary 8. Any adjoint singular reduction operators Q = 82 + C(^) u)8u and Q* = 81 + 
C,*{x,u)8u of equation ([5]) (where necessarily 7^ and C* 7^ are associated with the same 
one-parametric family of solutions of this equation. 

Let C be an arbitrary solution of equation Q. Then C„ 7^ and Q = 82 + C{x,u)8u is 
a reduction operator of equation ([5|). Consider a one-parametric family = {u = f{x,>c)} 
of solutions of ([5]), which are invariant with respect to Q. (Such a family exists in view of 
Theorem [5j) Expressing the parameter x from the equality u = f{x,x), we obtain that x = 
u), where 7^ 0. C = ^2 = for any u £ i.e., for any admissible values of (x, x). 

This implies that the representation C = —^2/^u is true for any admissible value of (x, u). This 
provides the background for the following statement. 

Corollary 9. The nonlinear three-dimensional equation ([6]) is reduced by composition of the 
Bdcklund transformation ( = —^2/^u, C* = —^i/^u, where ^ is a function of {x,u), and the 
hodograph transformation 

the new independent variables: xi = xi, X2 = X2, x = 
the new dependent variable: u = u 

to the equation ([5]) for the function u = u{x, x) with x playing the role of a parameter. 

Proof. We take an arbitrary solution of equation ^ (the condition 7^ is implicitly 
assumed to be satisfied) and set C* = (-^ ~ Ci)/Cm- Iii view of the Frobenius theorem, the 
system $2 + C^u = 0, <I>i + C*^u = with respect to the function <I> = <I>(x,n) is compatible 
since its compatibility condition (i + C*Cu = C2 + CC coincides with ([6]) and hence is identically 
satisfied. We choose a nonconstant solution ^ of this system. Then 7^ 0, = —^2/^u 
and (^* = — $i/<I>„. After the hodograph transformation, the latter equations take the form 
— Ci^iu) and Uxi = C*(^!^)- This directly implies that for any value of x the function 
u = u{x,x) satisfies equation The parameter x is essential in u since 'U;^^ =!/<!>„ 7^ 0. 

It follows from the proof of Theorem [5] that the application of the inverse transformations to 
a one-parametric family of solutions of equation ([5]) results in a solution of equation ([6]). □ 
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Note 4. For any equation from class ([5]) with = 0, reduction operators of the form Q = 
82 + C{x,u)du, where / (resp. Q* = di + C*{x,u)du, where Cu 0) ^-Iso are bijectively 
associated with one-parametric famihes of its solutions, having the form {u = f(x,>c)} where 
fi^ / (resp. /2>c / 0). The one-parametric families with fi;^ = (resp. /2>< = 0) necessarily 
existing in this case correspond to ultra-singular reduction operators with Cu = (resp. = 0), 
and the correspondence is not one-to-one. 

The above investigation of singular reduction operators of nonlinear wave equations of the 
form ([5]) shows that for these equations the natural partition of the corresponding sets of reduc- 
tion operators is into triples of subsets singled out by the conditions 

1)^1 = 0; 2)^2^0; 3) C'^V 0. 

After the factorization with respect to the equivalence relation of vector fields, we obtain three 
subsets of reduction operators, which have to be investigated separately. The defining conditions 
for these subsets are, respectively, 

1) = 0, = 1; 2) e = 0, ^1 = 1; 3) e y^O, f = 1- 

Since any equation from class ([5]) admits the point symmetry permuting xi and X2, the second 
case is reduced to the first one and can be omitted. Finally we have two essentially different 
cases after factorization: the singular case = 0, = 1 and the regular case 7^ 0, = 1. 
The gauge = 1 is not uniquely possible in the regular case and may be varied for optimizing 
the further consideration of this case. 
Consider the other standard form 

un - U22 = F{u) (7) 

of nonlinear wave equations, obtained from ([5]) via the point transformation xi = xi — X2, 
X2 = xi + X2, u = u. Using this transformation, all the results derived for class ([5]) can easily be 
extended to class ([7]). Thus, any equation of the form ([7]) possesses two singular sets of reduction 
operators, singled out by the conditions = — and £} = and one regular set of reduction 
operators, associated with the condition 7^ it^^. The singular sets are mapped to each other 
by alternating the sign of X2 and hence one of them can be excluded from the consideration. 
After factorization with respect to the equivalence relation of vector fields, we have two cases 
for our further study: the singular case = = 1 and the regular case 7^ ±1, = 1. 
For nonlinear wave equations of the general form 

Mil - {G{U)U2)2 = F{u), 

where G{u) > 0, the natural partitions of the sets of reduction operators are determined by 
more complicated conditions depending on the parameter-function G. We will not discuss these 
equation here. We only remark that the singular sets of the corresponding reduction operators 
are associated with the conditions = VG^,^ and .^^ = — respectively. 

The above examples underline that the application of the conventional partition for factoriza- 
tion of sets of reduction operators often leads to the splitting of uniform cases and to combining 
essentially different ones. As a result, the derived systems of determining equations for the 
coefficients of reduction operators is far from optimal and difficult to investigate. Therefore, 
natural partitions based on taking into account the structure of singular families of reduction 
operators offers a decisive advantage. 
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7 Reduction operators and parametric families of solutions 



Proposition 4. Let Q be a reduction operator of an equation C Then the weak singularity 
co-order of Q for C equals the essential order of the corresponding reduced ordinary differential 
equation. 

Proof. We carry out a point transformation in such a way that in the new variables the opera- 
tor Q has the form Q = ■ (For convenience, for the new variables we use the same notations 
as for the old ones.) Then an ansatz constructed with Q is u = (p{io), where (p = (p{u)) is the 
new unknown function and lo = xi is the invariant independent variable. The manifold Q(r) is 
defined by the system Ua = 0, where a = (ai, 02)) a2 > ^, Q.i + a2 ^ r = ordL. 

Since Q G Q(>C), there exist diff'erential functions \ = \[ip\ and L = L[ip\ of an order not 
greater than r such that -Z^|u=(^(w) = (cf. [26]). The function A does not vanish and may 
depend on X2 as a parameter. The function L is assumed to be of minimal order f which may 
be attained up to the equivalence generated by nonvanishing multipliers. Then the reduced 
equation C: L = has essential order f. 

The condition wsco£ Q = k means that there exists a strictly kth order differential function 
L = L[u] and a nonvanishing differential function A = A[n] of an order not greater than r, which 
depend at most on x and derivatives of u with respect to xi, such that L\q^^^ = XL\Q^^y 

If f would be less than k, we could use Anew = Mu-^<p and Lnew = L\u-^<p in the definition of 
weak singularity and would arrive at the contradiction wsco£ Q ^ ordLnew = f < k. Therefore, 
f ^ k. (Here, "y ~^ z" means that the value y should be substituted instead of the value z.) 

Suppose that f > k. We have the equality XL = (AL)|„=^(j^') in which the variable X2 plays 
the role of a parameter. Fixing a value of X2, we obtain the representation 



L = A[ip] L 



^ Mu=ip{uj) 



U=ip{Lij), X2=X2 



X2- 



Since ordLj^^^j-^-) .^.^^^o ^ k < f, this representation contradicts the condition that f is the 

essential order of the reduced equation £. Therefore, f = k. The inverse change of variables 
preserves the claimed property. □ 

Corollary 10. Let Q be a reduction operator of an equation C Then the weak singularity co- 
order of Q for C equals the maximal number of essential parameters in families of Q -invariant 
solutions of C 

Proof. The essential order f of the reduced ordinary differential equation C associated with Q 
coincides with the weak singularity co-order of Q for C The maximal number of essential 
parameters in solutions of C equals the order of C. The substitution of these solutions into 
the corresponding ansatz leads to parametric families of Q-invariant solutions of C, and all Q- 
invariant solutions of C are obtained in this way. Therefore, the maximal number of essential 
parameters in families of Q-invariant solutions of £ equals f. □ 

Corollary 11. Let Q be a kth co-order weakly singular reduction operator of an equation C. 
Suppose additionally that a differential function of minimal order, associated with L on the 
manifold Q(j,) up to a nonvanishing multiplier, is of maximal rank in the derivative of u of the 
highest order k appearing in this differential function. Then C possesses a k-parametric family 
of Q-invariant solutions, and any Q-invariant solution of C belongs to this family. 

Proof. Under this assumption, the reduced ordinary differential equation C associated with Q 
can be written in normal form and hence has a A;-parametric general solution which contains all 
solutions of C Substituting it into the corresponding ansatz, this solution gives a A;-parametric 
family of Q-invariant solutions of C. There are no other Q-invariant solutions of £. □ 
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Corollary 12. Suppose that a differential function of minimal order, associated with L on the 
manifold Q(,.) up to a nonvanishing multiplier, is of maximal rank in the highest order derivative 
of u appearing in this differential function. If the maximal number of essential parameters in 
families of Q-invariant solutions of C is not less than the weak singularity co-order of Q for C 
then Q is a reduction operator of C. 

Proof. Point transformations of the variables do not change the claimed property. We use the 
variables and notations from the proof of Proposition HI Consider the differential function 
L[ip] = i|M=<^(a;)- It depends on X2 as a parameter and ordL = k. Due to the condition of 
maximal rank, we can resolve the equation L = with respect to the highest order derivative 
ipik). ^pik) _ Ji^ipj^ where ordi? < k. 

If 7^ 0, splitting with respect to X2 in the equation L = results in an ordinary differential 
equation = of an order lower than k. Any Q-invariant solution of £ has the form u = (p(co), 
where the function satisfies, in particular, the equation R[ip] = 0. This contradicts the 
condition that the maximal number of essential parameters in families of Q-invariant solutions 
of C is not less than k. 

Therefore, = 0, i.e., the equation ip^''^ = R[f] is a reduced equation which is obtained 
from C by the substitution of the ansatz u = ^{to) constructed with the operator Q = 82- □ 

Note 5. For any operator Q, the maximal number of essential parameters in families of Q- 
invariant solutions of C cannot be greater than wsco^ Q. 

Summing up the above consideration, we can formulate the following statement. 

Proposition 5. Suppose that a differential function of minimal order, associated with the dif- 
ferential function L[u] on the manifold Q(^) (r = oidL) up to a nonvanishing multiplier, is of 
maximal rank in the highest order derivative of u appearing in this differential function. Then 
any two of the following properties imply the third one. 

1) Q is a reduction operator of the equation C: L = 0. 

2) The weak singularity co-order of Q for C equals k (0 ^k ^r). 

3) The equation C possesses a k-parametric family of Q-invariant solutions, and any Q- 
invariant solution of C belongs to this family. 

The properties of ultra-singular vector fields as reduction operators are obvious. 

Proposition 6. 1) Any ultra-singular vector field Q of a differential equation C is a reduction 
operator of this equation. An ansatz constructed with Q reduces C to the identity. Therefore, 
the family of Q-invariant solutions of C is parameterized by an arbitrary function of a single 
Q-invariant variable. 

2) If the family of Q-invariant solutions of C is parameterized by an arbitrary function of a 
single Q-invariant variable then Q is an ultra- singular vector field for C 

8 Reduction operators of singularity co-order 1 

Encouraged by the above investigation of evolution and, especially, wave equations, we study co- 
order one singular reduction operators of general partial differential equations in two independent 
and one dependent variables. 

Consider an equation C: L = 0, where L = L[v\ is a differential function of order r > 1. 
Suppose that the function L admits a first co-order singular module of vector fields. (In view 
of Corollary [H we can restrict ourselves to considering only strong singularity of vector fields 
for differential equations.) Without loss of generality, up to changing variables we can assume 
that the module contains a first co-order singular set S = {Q''} of vector fields in reduced form. 
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i.e., = ^di + ^2 + C^u for any smooth function C of {x,u) and a fixed smooth function ^. 
Additionally, we can assume ^ E {0, u}. 

By Theorem [H the differential function L can be written in the form L = L{x, ^r,i): where 

Qr,! = K = + D2r'u, Qi ^ 1, ai + as ^ r) , 

and 7^ for some Wq, with ai = 1. Then the restriction of L to Q^^^ coincides with the 
restriction, to the same manifold Q^^^j of the function L'^ = L{x, 0.r,i), where 

nr,i = (D°i(Q^)°^u,ai ^ l,ai +02 ^ r). 

Thus, the form of L'' is determined by the forms of L and ^ and a chosen value of the parameter- 
function Depending on the value of the differential function L'^ may either identically vanish 
or be of order or 1. This means that either the vector field Q'^ is ultra-singular or sco^ Q'^ = 
or scolQ^ = 1, respectively. We investigate each of the above cases separately. Below we 
additionally suppose that the function L'^ is of maximal rank with respect to u (resp. ui) if 
scol Q'^ = (resp. sco^ = 1). 

The values of C for which Q'^ for C is ultra-singular are singled out by the condition = 0, 
where u and ui are considered as independent variables. Splitting this condition with respect 
to ui gives a system 5_i of partial differential equations in C of orders less than r, which may 
be incompatible in the general case. The incompatibility of this system means that the set S 
contains no ultra-singular vector fields. For example, evolution equations of orders greater than 
1 and nonlinear wave equations of the form ([5]) with F„ 7^ 0, in contrast to equations of the 
form dS]) with Fu = 0, have no ultra-singular vector fields, see Sections [5] and [H C, satisfying the 
ultra-singularity condition guarantees that Q'^ G Qi'^) ^^'^ the family of Q*" -invariant solutions 
of C is parameterized by an arbitrary function of a single (^''-invariant variable. 

If scol Q'' = 0, the parameter-function ( satisfies the condition Z^^ = with u and ui viewed 
as independent variables, which is weaker than the ultra-singularity condition. Therefore, the 
corresponding system So of partial differential equations in of orders less than r, obtained by 
splitting the zero co-order singularity condition with respect to ui, has more chances of being 
compatible than Thus, any nonlinear wave equation of the form ([5]) with Fu ^ admits 

zeroth co-order singular vector fields although this is not the case for ultra-singular vector fields. 
At the same time, evolution equations do not possess zeroth co-order singular vector fields. 

Certain conditions which are sufficient for the compatibility of Sq can be formulated. Thus, 
if L^^^ „j = and = then the system Sq is compatible since it is satisfied by any ( with 
Cu = 0. In other words, sco^ ^ for any ( = Ci^)- Let us consider this particular case in 
more detail. (Recall that under the condition = the coefficient ^ can be assumed, up to 
point transformations, to equal but we will not use this possibility.) 

If additionally L^^^^ = 0, the condition L'^ = under the assumption ( = C,{x) implies only 
a single partial differential equation with respect to Any of its solutions is a solution of 5_i 
and hence the corresponding vector field Q'' is ultra-singular for L. 

Otherwise scol Q'^ = and we can resolve the equation L'^ = with respect to u: u = G'^(x), 
where the expression for the function depends on the parameter-function (" = ({x) and its 
derivatives up to order r — 1. Then the conditional invariance criterion is equivalent to the rth 
order partial differential equation = ^G^ + with respect to C- If C is a solution of this 
equation then is a reduction operator of C. The ansatz constructed with the operator 
can be taken in the form u = (p{oj) + G^(x), where ip = ip{u>) is the new unknown function and 
io = uj{x) is the invariant independent variable satisfying the equation + 1^2 = 0. It reduces 
the initial equation i2 to a trivial algebraic equation ip = 0, i.e., the function u = G'^{x) is a 
unique Q^-invariant solution of C. Conversely, let us fix a solution u = f{x) of the equation C 
and set C = C/i + /2. Then / = G<{x) and hence ( = ^G^ + G^, i.e., in view of the conditional 
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invariance criterion = ^di + ^2 + C^u is a reduction operator of C, and = 0. The solution 
u = f{x) is invariant with respect to by construction. Thus we obtain: 

Theorem 6. Suppose that an equation C: L = possesses a first co-order singular set S = {Q^} 
of vector fields in reduced form Q'' = ^c?i + ^2 + Cdu with ^„ = 0, i.e., its right hand side L is 
represented in the form L = L{x,ilr,i), where 

^r,l = (o^a = i^Di + D2)"'U, Oi ^ 1, Oi + 02 ^ r) , 

Lujc, / for some a with ai = 1, and additionally L^^^ o) ~ *-* ^'^'^ ^^(o a) ^^^^ there exists 

a one-to-one correspondence between solutions of C and reduction operators from S with = 0. 
Namely, any such operator is of singularity co- order and corresponds to the unique solution 
which is invariant with respect to this operator. The problems of solving the equation C and the 
exhaustive description of its reduction operators of the above form are completely equivalent. 

Now we consider the regular values of ( for which the singularity co-order of Q'^ coincides with 
the singularity co-order of the whole family S (and equals 1). If sco^^Q^ = 1, the parameter- 
function satisfies the regularity condition L^i 7^ 0. Therefore, the equation L'^ = which is 
equivalent to C on the manifold Q^^^ can be solved with respect to ui: ui = G^(x, u), where the 
expression for the function depends on the parameter-function (" and its derivatives up to 
order r — 1. Applied to the equation C and the operator Q^, the conditional invariance criterion 
implies only the equation 

Ci + CuC^ - iCi + e«G'^)G^ = CGi + Gi + CGi (8) 

with respect to the function C. Therefore, in this case the system of determining equations 
consists of the single equation ^ and, therefore, is not overdetermined. This equation can be 
rewritten as the compatibility condition 

Ci + CuC^ - ici + e.G'^)G^ - aoi + cic^) = + (c - ^G<)Gi 

of the equations ui = and ^ui + ^2 = C with respect to u. The order of ([8]) equals r and 
hence is greater than the order of the system Sq. This guarantees (under certain conditions 
of smoothness, e.g., in the analytical case) that the equation ([8]) has solutions which are not 
solutions of Sq. In other words, the equation C necessarily possesses first co-order singular 
reduction operators which belong to S. 

The results of Section [7] imply that for each first co-order singular reduction operator Q of the 
equation C there exist a one-parametric family of Q-invariant solutions of C. If the equation C 
admits a co-order one singular module of vector fields, the converse statement is true as well. 

Theorem 7. Suppose that an equation C: L = possesses a co-order one singular set S = {Q"^} 
of vector fields in reduced form Q'' = ^di + 52 + C^u- Then for any one-parametric family T 
of solutions of C there exist a value of the parameter-function ( = (^[x,u) such that Q'^ is a 
reduction operator of C and each solution from T is invariant with respect to . 

Proof. Consider a one-parametric family = {u = /(x, x)} of solutions of C. The derivative 
fj^ is nonzero since the parameter x is essential. From u = f{x, x) we derive x = $(x, u) with 
some function $ = ^{x, u), where ^ 0, and then define ^ = Ci^^ u) by the formula 

^ g^i + ^2 

Since fi = -{<^i/^u)\u=f, i = 1,2, then ^\u=ffi + /2 = CU=/, i-e-, any solution from T is 
Q^-invariant. Then either is an ultra-singular vector field for L or sco^ Q'^ = 1. (The case 
scoi Q'' = is impossible since otherwise the equation C could not have a one-parametric family 
of Q^-invariant solutions.) Any ultra-singular vector field for L is a reduction operator of C If 
scoi = 1 then Q is a reduction operator of £ in view of Corollary 1121 □ 
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Corollary 13. Suppose that an equation C possesses a first co-order singular set S = {Q^} of 
vector fields in reduced form Q'^ = ^di + ^2 + C,du, dn-d that no element of S is ultra- singular 
for C. Then up to the equivalence of solution families there exists a bijection between one- 
parametric families of solutions of C and its first co-order singular reduction operators belonging 
to S. Namely, each operator of this kind corresponds to the family of solutions which are in- 
variant under it. The problems of the construction of all one-parametric solution families of 
the equation C and the exhaustive description of its reduction operators of the above form are 
completely equivalent. 



This bijection is broken in the presence of ultra-singular vector fields. 

The above relation between one-parametric families of solutions and first co-order singu- 
lar reduction operators can be stated as a connection between the initial equation C and the 
determining equation ([8]). 

Corollary 14. Suppose that an equation C: L = possesses a first co-order singular set S = 
{Q'^} of vector fields in reduced form = £^di + ^2 + C^u- Then the determining equation 
for values of C corresponding to first co- order singular reduction operators of C is reduced by 
composition of the Bdcklund transformation 

+ $2 + C^n = 0, ^1 + G^^n = 

where ^ is a function of{x,u), and the hodograph transformation 

the new independent variables: xi = xi, X2 = X2, x = 
the new dependent variable: u = u 



to the initial equation C for the function u = u{x, x) with k playing the role of a parameter. 

Proof. We fix an arbitrary solution C, of equation ^ , which additionally satisfies the condition 
7^ 0. In view of the Frobenius theorem, the equations ^$i + <I>2 + C^?i = s-iid ^i-\-G'^^u = 
are compatible with respect to the function ^ = ^(x,u) since their compatibility condition 
coincides with ([8]) and hence is identically satisfied. We choose a nonconstant solution $ of both 
these equations. Then 7^ and 

After the hodograph transformation, the latter equations take the form ^Uxi +Ux2 = C(^) u) and 



u 



Xl 



G''{x,u). This directly implies that for any value of x the function n = n(x, x) satisfies 
the equation C. The parameter x is essential in tt since = l/^u 7^ 0. 

It follows from the proof of Theorem [7] that the application of the inverse transformations to 
a one-parametric family of solutions of the initial equation C results in a solution of equation ^ 
if the defined value of ( satisfies the regularity condition / 0. □ 
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